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ABSTRACT 


It is desirable to perform a study of the elastic 
responses of submarine hulls and internally attached 
structures to underwater shock loading. In view of the 
expense and possible involvement of non-shock testable 
equipment it is not always desirable to perform actual shock 
tests utilizing real hulls or SSTV's (Submerged Shock Test 
Vehicles). With the advent of large scale computing power, 
numerical methods now exist to predict equipment/hull 
responses. The ELSHOK, (ELASTIC SHOCK) computer code is 
used to perform a parametric study on a submerged shell with 
ameinternally attached substructure. Of particular interest 
is the phenomenon of dynamic amplification or lack of such 
mene responce or this Darticular substructure and effects 
on substructure-shell interaction due to varying | 
substructure mass and stiffness subjected to underwater 


explosion side-on shock loading. 
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l. See RODUCTION 


A. BACKGROUND 

A submarine must be able to withstand moderate to severe 
shock loadings that will result from underwater explosions 
(UNDEX) it aii ge be subjected to in modern warfare. 
Conventional and nuclear weaponry deliver devastating forces 
and no vessel can be expected to survive in the events of 
near or direct hits. The emphasis on design and testing 1s 
in the area of survival of mission critical equipment, 
machinery and weapons systems of a platform. The ability of 
a naval vessel to carry out itsS mission after being 
subjected to an UNDEX will depend on the survivability of 
these installations. 

The current specifications for building ships and 
submarines contain the requirements for shock loading to be 
met by the builder or vendor of shipboard installations 
and equipment. In general, all critical equipment is 
required to pass a series of shock tests where testing is 
practical as outlined in MIL-S-901D [Pef. 1}. This document 
specifies the shock testing requirements for shipboard 
machinery, equipment and systems which must resist High 
Impact (HI) mechanical shock. The purpose of these tests is 


to determine the suitability of machinery, equipment and 


systems for use after exposure to severe shock which may be 
incurred in wartime. It is not always practical -tto test 
equipment due to size and weight. In the case of equipment 
design, it is not always practical to evaluate an 
installation prior to final design. Mechanical shock 
testing of individual equipments may not always indicate how 
Tet ene Will perform or respond when actually installed. 
An analysis for a given design is specified for the 
situations where testing cannot be conducted. It is based 
on limited and dated testing and observations which are used 
to specify a shock desigin factor. A static design force is 
arrived at for various weights of equipment which is used to 
evaluate the design of support structures of installatiomes 
The method is simple to apply but does not always correlate 
well with actual testing. 

Numerical methods have been developed for use in 
submarine shock response analysis. The computer aa ELSHOK 
(Elastic SHock) developed by Weidlinger Associates for the 
Office of Naval Research [Ref. 2] is used me conduct this 


study. 





Bae PURE OSE 





ELSHOK has been developed through a series of controlled 


tests and has been used to predict equipment responses prior 


to low level explosive tests on submarines. Although ELSHOK 


has been validated by extensive testing it is considered 


10 





prudent to conduct additional research to gain insight to 
the sensitivity of the code's predictions when a variety of 
conditions are applied. The intent of this investigation is 
to examine the change in the equipment/hull Teraecions 
under these conditions. The models used to conduct the 
parametric study are Similar to those used in the mid 1970's 
in tests conducted by the Office of Naval Research (ONR) 
[Ref. 3]. The internal structure used was specifically 
designed to reduce the effects of dynamic amplification for 
the conditions applied in these tests. The same model 
geometry was studied here by varying one main dimension to 
tune the substructure to the natural frequencies of hull to 
see if this particular internal structure configuration 
retained this design feature over a range of mass and 


Stmpeness variations. 


et 


TI. UNDEX SHOCK RESPONSE USING ELSHOK 


A. ELSHOK PRINCIPLES OF OPERATION 

The ELSHOK computer code explained in [Ref. 2] 185 a 
Suite of computer programs assembled and developed to 
calculate the submerged shock response of a linearly 
elastic structure subjected to an underwater shock wave. 
Specifically, a ring-stiffened shell of revolution of finite 
length, with or without internal structure, iS considered to 
be immersed, initially at rest, in an@intinite acousme. 
fluid and to be excited by an acoustic shock wave. The 
shock wave originates from an UNDEX arbitrarily located from 
the structure and 1s assumed to propagate as a spherical 
wave eae the medium. This assumption allows incident 
fluid particle velocities impinging on the wet surface of 
the immersed structure, “ne hull of a submarine, to be 
determined from the incident fluid pressures. FELSHOK 
contains two possible representations of incident pressure. 
One uses a known pressure-time never ia. This type 1S Of Mime 
kind that may be used to describe a tapered charge simu- 
lating a nuclear explosion. The other is the empirical 
decaying exponential put forth by Cole [Ref. 4] for spheri- 
cal charges of conventional explosives. If the shell 


contains internal equipment, a substructuring technique 


eZ 

















described by Ranlet [Ref. 5] treats the internal equipment 
response by coupling the free-free modes of the empty ring- 
stiffened shell and the fixed-base modes of the internal 
equipment through use of dynamic boundary conditions. 
Structure-fluid interaction is calculated using the Doubly 
Asymptotic Approximation (DAA) presented by Geers [Ref. 6]. 
Normal fluid displacement of the Structure-fluid interface 
is expanded in a series of surface expansion functions which 
are orthogonal over the entire wet surface of the submerged 
hull or shell model. The surface expansion functions lead 
to matrices in which the elements are determined by matching 
the limits of zero and infinite frequencies of the pressure- 
velocity aeniaeionehaiess Thus, exact solutions are obtained 
‘for the transient problem at early and late times and DAA 
results in a smooth transition between the two limits. 
Because DAA calculates the eGitects Of them lumide in the 
interaction problem through the series of surface expansion 
functions defined only at the wet surface, the fluid field 
is in effect uncoupled from the structural field. The fluid 
effects produced by DAA are simply applied as additional 
loading in the modal component structural analysis of the 
shell. 

The substructuring ee ieee separately determines the 
vibration modes of the components that comprise the combined 


submerged structure of shell and internal substructures. 


eS 


Equations of motion of the entire structure are obtained by 
combining interaction forces and moments and enforcing 
compatibility of deformation at component connection points. 
This precludes calculations involving the modes and natural 
peecueneee of a system of combined components and also 
avoids a requirement for a combined system stiffness matrix. 
The ELSHOK code also handles the case of no substructure in 
the shell using the same substructuring method. 

ELSHOK is particularly suited to the task of a 
parametric study due to its purely numerical nature and 
modular organization. Various internal equipment models can 
be combined repeatedly with the same shell model. Inputs of 
presoumes time profiles or decaying exponential types are 
mathematically applied to these numerical models and result 
in output in the form of velocity-time histories. This 
transient response analysis takes into account interactions 
of shell and subStructure and eliminates dependency on data 


derived from actual explosion testing. 


B. ORGANIZATION AND IMPLEMENTATION OF ELSHOK 
The implementation of ELSHOK is accomplished uSing a 
Suite of seven major programs. These major components are: 
1) BOSOR4 - Structural analyzer for shell [Ref. 7] 
2) ACESNID - virtual mass processor 
3) PIFLASH - shell-fluid processor 


4) SAP IV - structural analyzer for subStructure [Ref. 8] 
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o) bf Perea l = substructure processor 
6) USLOB - time integration processor 
7) PUSLOB - plotting processor 
Illustrated in Figure 1 are the general relationships 
between the main components and four ese s Preo wih rer 
ELSHOK is separated. In any analysis using ELSHOK the shell 
‘model and internal structure model must be available or 
constructed by the user. Development of the models used in 
this study is discussed later. Once models are arrived at 
the analysis can be carried out in the four phases by 
executing each of the seven major programs sequentially in 
the order listed above. 
ie nee oral 1 and Fluid Anabys Us 

The finite element difference code BOSOR4 is used to 
analyze the idealized model of the submerged hull or shell. 
The model must be symmetric about its longitudinal axis and 
it may contain rings and bulkheads. Circumferential rings 
can be treated as discrete stiffeners or, if closely spaced, 
they may be represented in the orthotropic approximation 
(smeared) [Ref. 9]. An entire main body or full model is 
used to capture the gross effects of the shell response when 
calculations are made for circumferential wave numbers N = 0 
and 1. These include rigid body translation and eT opi ae on 
the structure (N = 1) and the torsional modes of the struc- 


Paw eter tOrming the analysis for N = 0 and > 2 


15 


(breathing mode N = 0) obtains the local shell response. In 
this analysis wave numbers N = 0, 1, 2, and 3 are used. 
Figure 2 shows the basic nodal patterns for N = 1-4. ELSHOK 
requires separate runs of BOSOR4 for each wave N used in the 
analysis. Each BOSOR4 run eerie the in-vacuo free-free 
modes and corresponding natural frequencies from the dynamic 
properties and geometry supplied in the model of the shell. 
If the shell is symmetric about the midpoint of its length 
and the shock wave impinges on it so as to produce symmetric 
reponse, only half of the structure need be modeled. Such 
1s the case in this study since the shock wave emanates from 
an Origin located normal to the shell at midships. 

| ACESNID (Accession to Inertia and Damping) accom- 
plishes the second step of Phase I by calculating the 
virtual mass array. This provides the late~time contri- 
bution of the Doubly Assymptotic Approximation [Ref. 6]. 
Normal displacements corresponding to the surface expansion 
functions are applied to the surface of a cavity in oHe 
fluid (infinite medium) having the same shape and size as 
the wet surface of the shell model being investigated. The 
cavity geometry 1s obtained from the shell geometry 
contained in the BOSOR4 output files. Only one run of 
ACESNID need be performed to encompass all values of N being 
used and all values of MN being used in BOSOR4 must be 


indicated in ACESNID for compatiblity. 


156 











PIFLASH (Prepare Input for Fluid and Shell) 
accomplishes the last step of Phase I and that is to 
reorganize the fluid and shell data files from BOSOR4 and 
ACESNID into a fluid-shell file that is used to solve the 
response equations. If there are no internal substructures 
involved in the analysis, then this fluid-shell file is the 
main input for Phase III and Phase II 1S not necessary. 


2. Phase II--Structural Analysis of Internal 
Substructures 


A finite element model for each internal 
substructure must be developed. The SAP IV code performs 
the structural modal analysis on each substructure to be 
included. The code contains a variety of types of elements 
that can be used with ELSHOK. Concentrated masses attached 
to the substructure may also be represented by the methods 
used in SAP IV. Each substructure requires a separate 
eco cence of calculations in which the fixed-base mode | 
Shapes and natural frequencies for each are determined. 
This data is recorded in the “substructure mode file" along 
with the geometry, stiffness and mass information of the 
model for future processing. 

PICRUST (Prepare Input for Calculating Response of 
PUMPsStLtuUcCtuLe) rearranges the data from the “substructure 
mode file" and computes additional information required to 
solve the response equations. PICRUST also provides for the 


selection of which modes are to be included from SAP IV and 


ey 


the specifics of connectivity between the shell and 
Substructure used in the response calculations. 
Additionally, influence coefficients corresponding to the 
base motions of the substructure, constraint modes and modal 
coefficients for calculation of interaction forces with 
moments between the shell and substructure are determined 
(Ref. 9]. PICRUST is the last step in Phase III and its 
results are stored in the substructure file. 

3. Phase III-~Shock Response and Time Integration 

The information required to complete the fluid- 

structure interaction is supplied as one input file. 
Therefore, if a substructure is being analyzed in the 
problem the shell-fluid file from Phase I must be combines 
with the substructure files produced in Phase II and become 
the input file for USLOB (Underwater Shock Loading of 
Bodies) which is ELSHOK's time integration processor. For 
the case of an empty shell the shell-fluid file is the input 
for USLOB. Integration is accomplished using a modified 
Runge-Kutta scheme. The time step and shock wave loading is 
specified by the user in USLOB. Care must be taken to be 
Sure that the product of the time step chosen and the 
highest circular structural frequency » of the modes being 
used from SAP IV is less than unity. Otherwise, the 
numerical integration scheme will not converge without 


changing the time step. 
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The cases studied in this work used both methods 
available to represent the incident fluid pressure. A 
pressure-time history profile shown in Figure 3 was supplied 
as a typical tapered charge simulation. Nine discrete 
pressure-time points are entered as input to describe the 
incident shock wave loading. Tapered charges are used in 
UNDEX testing to produce the same type of loading resulting 
from an underwter nuclear detonation. Here the impulse of 
the shock wave is delivered to the target over a longer 
period of time compared to the interval over which the 
impact of a conventional explosive charge is delivered. 

Both are assumed to emanate from an origin far enough away 
SO as not to include the erget in the gas bubble effects. 
ELSHOK is used for elastic response analysis only and hulls 
within the immediate vicinity of the gas bubble formed in an 
UNDEX would most likely have responses well into the plastic 
response range. Figure 4 shows the pressure-time profile 
used to represent the incident shock loading of a spherical 
charge of HBX-l, a SGeeare (omen explosive, of a weight which 
would deliver the same impulse over the same period of time 
as the duration of the tapered charge loading and detonating 
at the same range from the closest point of the target or 
standoff. This profile is an empirical exponentially 
decaying pressure-time history where the initial pressure 1s 


the peak incident pressure on the shell due to the charge 


Lo 


weight and standoff. It is seen from Figure 4 that the 
majority of impulse of the HBX-1 blast is delivered inside 
the first millisecond. The empirical equation over which 


USLOB performs integration to obtain the fluid particle 


velocities is of the form [Ref. 4]: 


K 


P,(Ryt) = K,(W/?/R) * exp(-t/8) (1) 
where, 
P,(R,t) = the incident pressure on the hull at range R (ft) 
from explosive source to point of interest (psi) 
Ky = multiplicative constant for incident pressure 
K5 = spatial decay constant for incident pressure 
ts = time after arrival of shock wave at point of 
interest (msec) 
W = weight of spherical charge (l1b.) 
~ K : 
6 = K wt73 (wl/3 pp) : time constant of : (2) 
O 5 . 
exponential decay (msec) 
K, = multiplicative constant for time constant of 
exponential decay 
Ky, = spatial decay constant for time constant of 


exponential decay 
USLOB creates a generalized velocity file for the 
shell-substructure and shell-substructure transformation 
file to be used in the last phase of ELSHOK in producamg 
plotted and printed velocity-time profiles. The trans- 
formation file contains information required to describe the 


motion at the shell-substructure connection points. 


Ze 





Ae eee == Plotting Velocity Profiles 

PUSLOB (Plots for USLOB) processes the results of 
the previous three phases so they may be plotted or printed 
as velocity-time histories. Velocity time-profiles are 
requested for various locations on the shell and 
substructure by specifying a call out point in the input 
files. The call out points are the nodes that are used to 
describe the models of the shell and substructure in BOSOR4 
and SAP IV. Selection of nodes is important because 
velocities can only be obtained at those locations. The 
shell velocities may be represented in either a global or 
local coordinate.system. The global 1S represented by X, Y, 
Z that locate the positive sense of X toward the bow or 
forward section of the shell, the positive sense of-Y is to 
port and positive Z is up. The local system is described by 
u-longitudinal (fwd pos) v-circumferential; w-inward normal 
from the shell. The substructure is always described by the 
coordinate system specified in its SAP IV finite element 


model. 
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Figure 2. Modal Patterns for Cylindrical Shells 
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III. MODELS INVOLVED IN THE ANALYSIS 


A. MODELING WITH FINITE ELEMENT/FINITE DIFFERENCE METHODS 
Finite element/finite difference methods have become 
powerful tools in engineering with the advent of large scale 
computing power. These methods are of great use in modeling 
techniques. Even though the models are much simplified for 
use in these techniques the results arrived at can be quite 
accurate or at least adequate for the problems that might be 

impossible to solve by other methods. Accuracy of the 
results is closely tied to how well the model represents the 
actual Reng eee oc, so care must be taken when 
constructing the model. In ELSHOK, the finite difference 
code BOSOR4 is used to describe the dynamic behavior of 
underwater hulls by representing them as much simplified 
models of ring-stiffened shells of revolution. Similarly, 
the models used to analyze internal substructures are also 
Simplified representations made with the finite element code 


SAP IV. 


B. SUBMERGED SHELL MODELING 

Figure 5 1S a schematic of the general shell structure 
and substructure model depicting the notation and convention 
used to describe the shell global coordinate frame X,Y,Z, 


shell local coordinate frame u,v,w and the substructure 


ZS 


coordinate frame x,y,z. The z-axis orientation of the logem 
Substructure coordinate system is related to the Z-axis of 


the shell by the angle a_ where o represents the substruc- 


0 
ture. Points are located along the meridian of the shell 
uSing the variable s. The angle ® is used to locate 
circumferentially the desired local shell displacements. 

One shell model is uSed in all calculations of this study. 
Once the symmetric full model was constructed no changes 
were required to conduct the investigation of varied 
Substructure sizes and types of UNDEX. 

The shell or hull that has been modeled for analysis by 
ELSHOK in this study 1S Similar to an actual structure 
COnencmaee for testing by the Office of Naval Research 
[Ref. 3]. Figure 6 is a schematic of the representation of 
the full shell model used. Though its size is much smaller 
it 1s also of Similar geometry to that of SSTV's (Submerged 
Shock Test Vehicles) that are used to investigate the 
underwater shock effects on submarine installed equipment. 
The model is an idealized submarine hull in that both can be 
conSidered as free-free ring-stiffened cylinders with end 
closures. The shell in this Study 1s in fact a cylinder of 
revolution made of high strength steel stiffened by six 
large steel rings plus forty small stringers also of steel. 
The end plates are aluminum. An internal substructure, also 


Similar to one involved in the testing by ONR is mounted at 
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midships on the port side of the hull. In BOSOR4 the model 
is divided into segments of consistent properties contained 
in each cross section. When the structure is symmetric 
Spout the midpoint of 1ts longitudinal axis only one half of 
the model need be configured. This is accomplished by 
specifying the appropriate boundary condition at the point 
of symmetry. In this problem symmetry is employed allowing 
the shell model to be divided into two segments. Segment 
one is the aft aluminum end plate and has nine nodes 
specified on its surface. Node 1/01, node one on segment 
one, is the center of the end plate. The other nodes move 
as concentric circles to the outer circumference. Segment 
two is the cylinder from the end ola to the midships 
point. Twenty-seven nodes are specified for segment two 
that are actually bands around its circumference. Node 
numbers increase from aft to forward such that node 2/27 is 
the midpoint on the cylinder. | 

Concentrated masses and substructures internal to the 
shell must be modeled separately if they contribute signi- 
ficantly to the stiffening of the shell. Therefore, the 
large rings have been modeled as discrete stiffeners and the 
internal substructure is modeled separately using SAP IV. 
If there were no interest in the substructure response it 
could be modeled as a concentrated mass attached to the 


shell. There is only one substructure,a doubly overhung 


ZT 


cantilever beam, used in this study and there are no 
concentrated masses. The smaller stringers do not 
contribute as significantly to the response of the shell as 
do the large rings so they are smeared onto the shell to 
account for their mass [Ref. 9]. 

The cylindrical section of the shell is a surface of 
revolution having consistent properties throughout its 
entire length. This allows information regarding Young's 
modulus, Poisson's ratio, mass density, thermal expansion 
coefficients, plating thickness and whether or not rings are 
Smeared to be specified at only one point for the segment in 
BOSOR4. Similar information is specified for segment one, 
the end plate. Tables I and II are listings of the natural 
Frequencies and wave numbers for each mode calculated. 
Twenty-two N = Q breathing modes, thirty N = 1 rigid body 
and whipping modes and thirty N = 2,3 modes were retained 
for analysis. The torsional modes for N = 0 were not 
retained to eliminate a rotational input at the base of the 
Substructure model. Modal contributions corresponding to N 
= 2,3 normal, meridional and circumferential displacements 
were retained for this full model case. Figure 2 shows the 


basic shapes of the waves on the cylindrical cross section. 


Ge INTERNAL SUBSTRUCTURE MODELING 
In ELSHOK calculations, a substructure model is formu- 


lated to determine the modes, masses and natural frequencies 
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of the internal equipment being analyzed. Phase II of 
ELSHOK utilizes the SAP IV finite element code to accomplish 
modeling of submarine-installed equipment. ELSHOK is 
written in program modules that execute independently and 
exchange data through the output files. This implementation 
lends itself well to the type of parametric study conducted 
in this thesis work. A variety of different internal models 
may be combined with the same shell model to gain insight on 
the interaction of the response between the shell and 
various substructure characteristics. Modeling the 
substructures in a finite element method more accurately 
predicts effects of the substructure on the hull to which it 
is attached than if the model is represented as a 
concentrated mass. A concentrated mass model would not 
account for any substructure responses that might become 
dynamically amplified. 

The substructure modeled for this study is a doubly 
overhung steel beam mounted with steel brackets to the port 
Side of the shell. Connection points are at the two large 
rings at frames 18 and 27, placing the beam symmetrically 
about the midships section. Only half of the beam structure 
requires modeling since symmetry may be employed. Symmetry 
1s indicated in SAP IV by specifying the appropriate 
boundary condition at the point of symmetry on the model. 


This 1s Similar to the procedure in BOSOR4. The dimensions 


Zo 


and details of the model are shown in Figure 7. Figure 7 
also shows the orientation of the substructure within the 
shell and to the shockwave. The doubly overhung beam is 
represented using the three-dimensional beam element model 
of SAP IV depicted in Figure 8 which shows half the 
structure with the first mode for D = 4.2 in. superimposed 
on it. In SAP IV nodes are used to divide the model into 
separate elements and serve the purpose of call out points 
for velocity-time histories in Phases III and IV of ELSHOK. 
Modeling one half of the doubly overhung beam was 
accomplished using 21 nodes, 4 on the bracket and 17 on the. 
beam half length. The brackets are represented as purely 
extensional members (a spring with mass) and the intro- 
duction of a bending moment at the connection between the 
beam and bracket is eliminated. This is accomplished by 
specifying an end release node at node 12 that represents 
the connection point on the eT The interest in the 
response of the substructure is limited to the z-direction 
and restrictions on the degrees of freedom of the model were 
appropriately specified. This reduces the calculations 
required and produces responses in the direction of 
interest. 

In order to investigate the response interaction between 
the shell and substructures of various weights and charac- 


teristics the depth D of the beam is changed for.each case. 


one 





It should be noted that the width of the bracket varies with 
the depth of the beam by 1.5 x D. Developing the various 
sized models was accomplished by editing the SAP IV input 
code to reflect the new cross-sectional areas and moments of 
inertia for the beam and bracket elements. Other PO Pei emo om 


remained unchanged. 


D. SHELL/SUBSTRUCTURE RESPONSE MODELING 

The structural problem of ELSHOK is separated into two 
parts to avoid calculations involving the modes, natural 
frequencies and stiffness of a combined system of the shell 
and substructure. Even though separate modal analysis 1s 
carried out on each model in different phases of ELSHOK the 
responses determined for the shell and substructure are 
arrived at by coupling the free-free modes of the shell with 
the fixed-base natural modes of the substructure. In effect 
the shell and substructure can be approximately represented 
by a two-degrees of freedom mass-spring model as in Figure 
9. Here the mass of the shell is M2 and the stiffness 
characteristics are K2. The substructure model is 
represented by Ml and Kl. It has been shown for some 
Situations uSing two single-degree of freedom models can 
produce nearly the same responses aS a two-degrees of 
freedom model when the mass ratio is very small. This in 
general does not hold for systems of large mass ratios and 


1S apt to produce marked differences in response between the 


Si 


two methods if the same stiffness characteristics are used 

in both [Ref. 10]. ELSHOK results can be viewed with more 
confidence in predicting component responses because it | 
utilizes the multi-degrees of freedom type model. It is 


preferred in this study presented due to the large mass 


ratios that are involved. 
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IV. ANALYSIS 


A. OBJECTIVES OF ANALYSIS 

The objectives of the analysis procedure are to 
investigate effects on component responses brought about by 
varying the mass and stiffness of the internal substructure. 
More specifically, it is divided into two major and one 
minor task. First of these is to determine how heavy the 
internal structure must be before it significantly affects 
the motion of the shell. Secondly, the dynamic amplifica- 
tion of the responses of the substructure and its effects on 
the sien are investigated. Lastly, it is of interest to 
know the magnitude of the forces generated in the support 
brackets. An integrating scheme is used to determine 
deflections and therebyses timnarve = ne teece ss 

The shell and substructure are subjected to side-on 
shock loading from an UNDEX. The UNDEX is represented in 
two different manners. One iS an incident pressure-time 
history simulating a tapered charge or nuclear detonation 
and the other is an empirical exponentially decaying 
incident pressure Simulating a spherical conventional charge 
detonation. The standoff or range from the origin of the 
UNDEX to the closest point of the target is the same for 


both the tapered and conventional cases. The weight of the 


conventional charge is selected to deliver the same impulse 
Ms that of the tapered charge. First, the transient 
responses of the empty shell are determined for a basis of 
comparison. Next, responses for a variety of internal 
Structure sizes are calculated for both conventional and 
tapered charge cases. The results are produced in the form 
of velocity-time histories and are obtained for points 


Specified on the shell and substructure. 


Be SONALYSIS PROCEDURE 

Geeescnel] natural Erequencies and modes for wave numbers 
N = 0, 1, 2, and 3 are found using the BOSOR4 code. The 
virtual mass array computed in ACESNID is done including 
each wave N specified in BOSOR4 for compatibility. In order 
to Simplify the analysis and since the brackets are designed 
to resist movement out of the horizontal plane, the 
torsional and rolling modes are eliminated from considera- 
tion to prevent a rotational input to the base of the 
Substructure. This is accomplished by setting the variable 
NTORSN = 0 and omitting any purely rolling or torsional 
modes for N = 0 from input in the PIFLASH code. All thirty 
modes for the other values of N are used. This shell model 
is used throughout the study in combination with each 
substructure. The BOSOR4 modes deleted and retained are 
indicated in Tables I and II. 

SAP IV was used to determine the fixed-base natural 
modes and frequencies as D was varied on the substructure 


a 


model. The investigation begins with a beam depth D of one 
inch and increases to twenty inches. It should be noted 
that twenty inches is nearly physically impossible to fit 
within the confines of the shell and was computed to 
investigate the effects of a large mass on the shell. Each 
SAP IV run produced ten fixed-base natural frequencies 
corresponding to different mode shapes of the doubly over- 
hung beam. Table III contains frequencies computed for the 
cases D = 1.0, 4.2 and 9.0 in. Not all the modes determined 
in each case were used in the follow-on calculations because 
some of the high frequencies would not allow convergence of 
the integrating scheme used without changing the time step. 
All modes having frequencies of 3979 Hz and greater were not 
used in the analysis. These higher modes do not contribute 
significantly enough to warrant varying the time step to 
remain within the convergence limits. 

Results from each SAP IV run require processing by the 
PICRUST code for preparation to merge the shell and 
Substructure data files with the fluid data files. The 
configuration of each substructure installation along with 
the modes to be used in Phase III calculations are 
specified. The substructure z-axis is related to the shell 
coordinate system by the angle O shown in Figure 5. The 
meridional locations of the attachment points are also 


specified here. Each case investigated placed the 


4() 
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Substructure attached to the port side or an angle of 180° 
from the shell Z-axis orientation (exactly one half the 
distance between the top and bottom). The brackets connect 
to the rings at frames 18 and 27 placing the substructure 
evenly about the midships section of the shell or node 2/27. 
Each case D is subjected to two types of UNDEX representa- 
mrons Naving Origins at 70 ft. (840 in.) from the shell 
abeam the starboard side at midships. The origin also is 
located in the shell's Z-X plane. After PICRUST prepares 
the shell/substructure-fluid data the time integration of 
the incident pressures experienced by the shell surface is 
performed to determine velocities. 

The USLOB code models the type of UNDEX by differently 
representing the incident pressures. The tapered charge 
pressure-time history is represented by inputting discrete 
pressure-time points that describe it. The conventional 
Bonerical charge is represented by equation (1) and 
specifying the appropriate constants for the type of 
explosive used. The weight of the charge specifies the size 
of the explosion and its resultant impulse and incident 
pressure on the shell. In determining a charge weight that 
would deliver the same impulse as the tapered charge 
pressure-time history equations from Cole [{Ref. 4] and 
constants describing HBX-l were used. Describing the 
impulse of unit area of the shock wave front up to a time t 


after its arrival at the shell wet surface by: 


4l 


te 
{ Pi(Reee (39) 
0 


LT (9 


where: P(t) P exp(t/%) (4) 


O 


Then integrating equation (2) over the time after arrival of 


the shock wave results in: 
Tit) 9Pse (= exni=o7aam ey 


Subs®rtuting eq@uatron (27) sar ae and re where, 


K 
. OS ae 
Py = K, (we/"/R) (6) 


into equation (5) then setting ie equal to the area under 
the curve in Figure 3 the charge weight arrived at through 
iteration is 352 lbm of HBX-l. The constants for HBX-1l 
are determined from data recorded in feet and msec. This 
required Ky and Ko to be adjusted to represent the empirical 
relationships in inches and seconds in order to be con- 
Sistent with the units used in this application of ELSHOK. 
The standoff and lateral placement of the origin of the 
UNDEX is Specified in USLOB along with the duration of 
integration and the time step to be used. A time duration 


5 


of 10 msec and a step of 3xl0 ~ sec was used to capture the 


peak response in each case. Also no effects due to surface 








cutoff or reflection of shock waves were considered. MThis 
Situation represents a deep submerged test simulation where 
bulk cavition due to the surface doesn't occur and there is 


no bottom bounce or surface cutoff effects. 


oo PLOTTING PROCEDURE 

Plotting of velocity profiles is accomplished by PUSLOB 
in the last phase of ELSHOK. Each profile is a record of 
the velocity history for a single node specified on the 
shell or substructure. The plot is limited to two angular 
locations of the shell nodes and velocities are expressed in 
in/sec due to the units used in this study. The version of 
ELSHOK ‘available at the Naval Postgraduate School is 
installed on the VAX/VMS system and plots were produced on 
Tektronix devices. In order to aid in comparing the 
velocities of different nodes the data files used to create 
the PUSLOB:plots were processed using program PLOTCONV. 
Program PLOTCONV rearranges the data files and converts the 
units of velocity from in/sec to ft/sec. This facilitates 
the use of plotting routines EASYPLOT or DISSPLA available 
on the school's IBM system and expresses the plots in the 
customary units of the shock community: velocity in ft/sec, 
time in msec. Program FISH] is a program developed for use 
with the DISSPLA plotting routine. Copies of programs 


PLOTCONV and FISH] are found in Appendix B. 


43 


D. FORCES IN THE SUBSTRUCTURE SUl rer. 

Since ELSHOK produces results in velocity-time 
histories, force calculations require some post processing. 
Forces generated in the support brackets would be mainly due 
to extension and compression because the bracket design 
resists motion out of the horizontal plane. Due to this and 
to simplify the analysis, inertia of the brackets is 
neglected in the force calculations, though their mass 1s 
included in the shock response calculations of ELSHOK. 
Program FORCES found in Appendix C was used to convert the 
velocity histories to deflections using a Simpson's one- 
third integrating scheme. Uniaxial forces estimated are 
based on F = ceEA where € is the estimated deflection of the 
bracket divided by the length between the connection points 
of bracket at the ring and beam junction. E is the modulus 
of elasticity of steel and A is the cross section area of 


the bracket, i... 5 See: 
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V. JRESUETS 


Numerous cases were explored in this parametric study in 
which the beam depth D varied over the range of one to 
twenty inches. ELSHOK generated a large amount of 
information about each case. The most pia ciesee results 
have been summarized in tabular and plotted forms which 
contain information on cases subjected to the tapered charge 
and the mentee charge UNDEX simulations. SAP IV 
predicted frequencies and modes for the tapered charge case 
with beam depth _— 4.0 ingaawee compared very closely with 
the ONR model previously tested. This provides a certain 
amount of confidence in the SAP IV model used throughout 


this investigation. 


A. SHELL AND SUBSTRUCTURE INTERACTION 

Even though the variation of substructure sizes included 
the cases resulting near the 600 Hz secondary frequency that 
was expected to induce dynamic amplification in the original 
ONR model there is no such response over the entire range of 
cases investigated. The case for D = 4.2 in. conventional 
charge resulted in a peak of velocities at the tip and 
midsection of the beam as shown in Figure 10 and Figure oe 
However, this 1S not a resonant condition. This slight 


amplification can be traced back to shell modes 1, 2, 3 for 


48 





N = 2 and mode 1 for N = 3. Tables IV and V summarize the 
range of results for both charge types and, when viewed with 
Figures 12 and 13, there is some indication of transfer of 
energy between the substructures and the hull, but not much. 
There seems to be no significant interaction between the 
beam and shell. The basic shape of the velocity-time 
histories at the beam tip, midsection and junction remain 
nearly the same throughout the range of D studied. The only 
effects seem to be the decrease of beam and bracket maximum 
velocities as the weight of the beam is increased. The 
complete velocity profiles for the beam bracket and connec- 
tion points for the D = 4.2 in. cases are shown in Figures 
14 through 17. This size was selected for comparison 
because it did demonstrate slightly higher overall velo- 
cities on the beam and also falls nearly into mid-range of 
the cases 1f D = 15.0 and 20.0 in. are regarded as unreason- 
able models due to their size. These two cases were looked 
at simply to check the effect of the large mass they would 
represent on the hull. Noting Figure 13 that shows the 
maximum velocities of the bracket connection points, the 
rings at frames 18 and 27, it is seen that the variation 

in peak responses is only slight over the range of D = 
4.0-10.0 in. for the tapered charge case. This corresponds 
to substructure to shell mass ratios of .07644-.19875. The 


empty shell velocity at this point is 10.48 ft/sec and its 
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highest maximum velocity occurs when D = 7.0 in. which 1s 
11.88 ft/sec corresponding to a mass ratio of .13638. The 
conventional charge empty shell velocity at this point was 
32.14 ft/sec and this was the maximum velocity recorded on 
frames 18 and 27 of the conventional cases. 

The velocity-time histories of the beam and bracket for 
the cases D = 2.0-10.0 in. for both charge types are very 
Similar to the D = 4.2 in. cases shown in Figures 14 through 
16. At D= 1.0 in. the response was very oscillatory and 
produced the largest maximum velocity for the conventional 
charge at the beam tip. This could be due to the large 
impulse combined with the smallest mass ratio of .01874. 

The other cases simply show overall decreases in velocities 
but retained the baSic shape of the velocity-time histories. 
Comparing the empty shell to the shell with Substructure 
velocity-time histories at the midsection of the shell and 
frame 18 demonstrates wiffit little Veemerr 1 bDutaoT Sere 
substructures' responses made in the shell response. In 
Figures 17 through 20 the shape of the velocity-time 
history profiles for these points can be seen to be mostly 
unchanged by the presence of the substructure. This was 
true for all cases, conventional and tapered, and increasing 
the weight of the substructure served mostly to reduce the 
velocity at the connection point. 

Comparing the frequencies of the shell model from Table 


II with the first five natural frequencies of the substruc- 
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ture models there are very few matched frequencies at and 
below 600 Hz which are the range expected to excite the 
first two modes of the substructure. The frequencies for 
the cases listed in Table III show how quickly they increase 
to magnitudes where their contributions to response is 
minimal. Any modes of the substructure being excited by the 
pase input at the connection points of the brackets and 
shell did not result in a dynamically amplified response. 
Each iecement of substructure size produces a set of fixed 
base frequencies of larger magnitude than the previous set. 


If the natural frequency is equal to: 


where K = stiffness and is = effective mass then the 
stiffness must be increasing more rapidly than the mass of 


each structure to produce the higher frequencies. 


B. OVERALL SHELL RESPONSES 

The shell response overall seems to exhibit some of the 
expected modes for the circumferential waves used in the 
calculations with the exception of the purely torsional 
modes that were deleted early on. First looking at the 
velocity-time histories in Figures 21 through 24 the extreme 
difference in velocities produced by the two types of 


charges 1S apparent. As expected, the conventional charge 


ol 


produces larger velocities in each case. These velocities 
are in the athwart ship direction, positive sense to port, 
at call out points on the midships station and the three 
frames where the discrete rings are located. The velocity- 
time histories due to the conventional charge illustrate how 
the velocities propagate from the midships station to the 
end plates. Looking at the peak velocities the shell can be 
seen to be bending about the midships station while the 
enero shell translates in the'direction of the shock wave. 
These modes can be seen in the velocity profiles for 
conventional and tapered charge cases where the midships 
peak velocity is less than those for the two rings located 
at frames 9 and 18. Looking at Figures 22 and 25 shows that 
the velocities at the midships station for the conventional 
case on the port and starboard sides are of opposite signs 
at about time equal to 1.0-2.0%msec. "fhis 1PPUStrarcom 
higher order circumferential wave of N = 2. Pilg ene and 
26 show the velocity profiles at the top of the shell in tne 
athwart ship direction. In each case run the top and bottom 
velocity profiles were exactly matched indicating the 


Symmetric response to the shock loading in the translational 





and whipping modes. Figures 27 and 28 illustrate the slight 
lag in the velocities on the port side. As a rule the 
conventional starboard side velocities at ring frames were | 


always higher than the port. This is due to the short 


SIZ 








Gumation Of shock front of the conventional Charge. 
Alternately, the tapered charge caused slightly higher peak 
velocities on the port side but also lagged the peak 
responses on the stbd side. This shows how the tapered 
charge allows the port side to exhibit more response to the 
eee@rting forces. 

The end plate velocities are shown in Figures 29 and 30. 
The end plate in Figure 29 can be seen to bulge in and out 
at the center for the conventional case and in Figure 30 the 
tapered case shows the same response. The ring at frame 
zero tends to hold the outer circumference of the end plate 
in place as shown by the lower velocities experienced at the 
outer edge. These velocity-time histories are indicative of 


breathing modes of the shell. 


C. FORCES IN THE SUPPORTS 

The axial forces were estimated for the D = 4.2 in. 
model subjected to both the conventional and tapered 
charges. In both cases the forces fluctuate between 
compressive and tensile as indicated by the changes in sign 
in the results in Appendix C. The maximum force estimated 
occurred in the conventional case. It was a compressive 
force of about 103 kips or a stress of about 16 psi. The 
maximum tensile force experienced by the support bracket for 
this case was estimated at approximately 75 kips. The 


tapered charge case produced much lower force estimates. 


DO 


The maximum was tensile and occurred at about 1.7 msec 
whereas the conventional maximum force occurred at about 

-3 msec. The maximum estimated tensile and compressive 
forces for the tapered charge case are 25 kips and 23 kips, 


respectively. 
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SUMMARY OF RESULTS--CONVENTIONAL CHARGES 


D(an) Pichi eer 2@Ekz >) MAXIMUM VELOCITIBS (FPS) M1/M2 


TIP(04) MID(21) JCN(12) 2/21(180°9) 


1 25.46 91.63 77.21 43 40am somo 3 ieeG .01874 
2 50.92 162.2) 55.66 41 .32mNN 37 599meD 20mm 0979 
3 76.38 240.6 51.59 39.90 33.44 27.85 .05696 
4 101.8 315.9 58.5] 48.48 30.77 28.02 .07644 
4.2 106.9 330.5 59.16 9 46. oqm2o .cQmo7ce .0765 
5 127.3 387.1 58.04 48.76 26.56 26.43 .09602 
6 152.8 453.4 54.42 47.06 24.18 23.45 .11615 
7 178.28 514.0 49. 20d RON 4) omen dee 13638 
8 203.7 568.7 43.33 38.3525 0lnmennocmar .15685 
9 228.7 616.3 37. 7e0—eecc. o1mNNs4 >> mmr se .1775 
10 254.6 659.6 32.35 29.65 23.48 22.90 .19875 
15 381.9 798.5 20.91: 20.34 19.65 19.69 . 30m 
20° 50902 “863.3 “ia. 17a ees .4219 


(Shell mass M2 = 3420.0 lbm, 


Empty shell v(2/21) = 34.14 FPS) 
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Figure 14. Velocity-Time Histories at Beam Tip for D = eZ Ln 
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Figure 15. Velocity-Time Histories at Beam Midsection for 
Die line | 
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Figure 16. Velocity-Time Histories at Beam/Bracket Junction 
fei) = 4,2 In. 
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Figure 18. Velocity-Time Histories at Frame 18 for Emoty 
Shell 
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Velocity-Time Histories at Frame 18 for 
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Figure 27. Velocity-Time Histories at Frame 18 of ube nny 
Shell from Conventional Charge 
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Beguce 26. Velocity—Time Histories at brame 18 of Empty 
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Figure 30. Velocity-Time Histories at End Plate for Tapered 
Charade 
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VI.  CONCEUSIGHS 


The results from the cases studied tend to show that a 
model can be designed to minimize dynamic amplification. 

In no case did a truly resonant situation occur. As the 
substructure model increased in size and mass so did the 
stiffness and the effect was as if a lumped mass was 
attached to the shell. Even though substructure sizes 
having frequencies that matched those of the hull were 
tested, the corresponding mode shapes being excited did not 
result in dynamic amplification of the beam responses. The 
beam model having a secondary frequency of 330.5 Hz became 
slightly amplified but only for the conventional explosive 
analysis. 

Though the models used in this study are based on real 
configurations it is important to keep in mind oe iene 
numerical study on a simplified model of an already simple 
structure. Even though the substructure design used did not 
become highly excited when combined with the shell in this 
study, the beam's behavior might change if combined with a 
Shell that has a frequency range and mode shapes that are 
different from those of the shell used in this study. The 
shell model was a relatively small structure that resulted 
in a majority of frequencies too high to excite the 


Substructure significantly. 
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ELSHOK has the capability of representing much more 
complex models than those used here. It has the capability 
of simulating a full sized submarine hull with internally 
attached equipment. The ELSHOK code has been shown to give 
fairly accurate results, especially in the early portion of 
shock Mieiens, through comparisons of its predictions to 
physical test results. However, like any numerical device 
it needs such comparisons to verify its results. On the 
other hand, conducting tests using ELSHOK has the distinct 
advantage of eliminating risk of damage to valuable assets. 
Once a model has been verified to satisfactorily predict the 
Prototype responses, a variety of internal equipment models 
could be applied to perform a variety of studies at 
fractions of the cost of building real models. Again the 
drawbacks are the pitfalls of any numerical endeavor, 
namely, the output is only as good as the input. Simplified 
models may a ins i fied results and poor assumptions 
must be guarded against. 

It 1S important to realize that the pasate GO fea s 
study were arrived at through the implementation of a 
potentially powerful tool for future research at the Naval 
Postgraduate School, namely ELSHOK. Further studies in the 
area of UNDEX using the code should be carried out to prove 
its validity or develop it into more useful applications of 


design and research. Suggestions on improvements to develop 


is 


the line of research in this thesis from a numerical 
sensitivity study into a realistic problem are to obtain or 
construct a representative full-size submarine model or 'seu7 
for BOSOR4 and model all internal equipments to be tested in 
appropriate detail while applying fewer restrictions to the 


degrees of freedom. 
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APPENDIX A 


INPUT FILES FOR ELSHOK 


1. BOSOR4 INPUT DATA FOR ALL CASES 

The following data file was used to generate the 
cylindrical shell model. Four separate runs of BOSOR4 were 
made corresponding to N = 0,1,2,3 by setting NOB, NMINB and 
NMAXB equal to each case N. All other inputs remain the 


same for each case. 


8 3 


RING-STIFFENED en aoe FLAT ENDPLATES--SYMMETRY USED--N=0 


2 NDIC = ANALYSIS TYPE INDICATOR 

1 NPRT = OUTPUT OPTIONS ( }aMeN IMUM = MEDIUM 3=MAXI MUM) 

0 ISTRES= OUTPUT CONTROL (O=RESUL ints =SIGMA, 2=EPSILON) 
z 2 EG = NUMBER OF SHELL SEGMENTS (LESS THAN 53? 
H SEGMENT 


NUMBER i i ji 
NMESH = NUMBER OF NODE P INTS $ (5.5 = MIN. : 98 me 1) 
3) FOR NODAL He T SPACING 


NTYPEH= CONTROL INTEGE ER 
HAPE= INDICATOR (1 oh ROR REOMESR OF MERID 
DIUS T° BEGINNIN OF SEGMENT (SEE P. 6 
ARTAL cOORB ESA MS AT BEGINNING OF SEGMENT 
US AT END OF SEGM 


ENT 
iE CTLC OF SE ee 


P= 5 
NTYPEZ= CONTROL (1 OR 3) FOR RE DCE oN AoE SOME tL rtoN 
ZVAL. = D FROM LEFTMOST SURE TO REFERENCE SURF. 
S ) Rees). ETC: R THIS SEGMENT? 
MAX= 20) OF D OGRE: E RINGS IN THIS SEGMENT 


Sreeeds 


NRINGS=_ NUMBER 
K=ELASTIC FOUND 
LINTYP= INDICATOR 
NLTY PE=CONTROL 


E. LB/IN**3)IN THIS SEG 
1, 2 OR 3) FOR TYPE OF LINE LOADS 
2,3) TYPE OF SURFACE LOADING 


NWALL= INDEX Tes ay, 5, 6, 7, 8) FOR WALL CONSTRUCTION 
0. 1080000E+08 8 E = YOUNG'S MODULUS’ FOR SKIN 
0: 3200000 U POLS SON'S RATIO FOR RIN f 
0. 2535000E-03 SM_ =MASS DENSITY: OF SKIN (E. a ALUM. =.00025 LB-SEC**2/IN**4) 


ANRS 
SUR CONTROL 
DO YOU WANT TO PRIN OUE 


MODDOOMOOOO ODWOOMOO-rH Ws) 
' 


COEEFELPELENT OF 
CONTROL ¢ OR hi 


RMAL_ EXPA NSION 
OR a FOR ADDITTON OF SMEARED oTIFFENERS 


yee 6 peut UT OR 1 Unk 


N THE JI AT T MERIDIONAL’ STATIONS? 
N DO YOU WANT TO PRINT OUT nay D LOADS ALONG MERIDIAN? 
H SEGMENT NUMBER 2 2 2 2 
25 NMESH = NUMBER OF NODE P ITS Ss 65 = MIN. : 98 Max, )( 2) 
3 NTYPEH= CONTROL INTEGER OR 3) 4--OR NODAL POLNT SPACING 
iL NSHAPE= INDICATOR (1,2 0 ae REOR GEOMETRY OF MERIDIAN 
16. 81250 Rl = RADIUS AT BEGINNING OF SEGMENT (SEE P. 6 
0 Za = AXIAL COORDINATE AT BEGINNING OF SEGMENT 
16. 81250 R2 = RADIUS AT END OF SE 
129. 0625 aD = AXIAL COORDINATE AT END OF SEGMENT 
0 IMP = INDICATOR FOR IMPERFECTION (9 = 
3 NTYPEZ= CONTROL (1 OR 3) FOR REFERENCE SURFACE LOC TION 
0. 1875000 = DISTANCE FROM LEFTMOST SUR BE TO REFERENCE SURF. 
N DO YOU WANT TO PRINT OUT of! S) {$) e tire SEGMENT? 
3 NRINGS= NUMBER MAX=20) OF SCRE E RINGS’ IN THIS SEGMENT 
2 NTYPE = CONTROL FOR IDEN TLFICATTON OF RING LOCATION (2=Z, 3=R) 
3. 500000 I) = AXIAL COORDINATE OF ITH RING, Z( l 
53.12500 Zit) = AXTAL GOORDINATONOM MEG RING? Zia 
103. 7500 I} = AXIAL COORDINATE OF ITH RING; Z( 3 
1 NIYPER= TYPE (0 OR 1 OR 2 OR 4 OR 5) OF DISCRETE RING NO.( 1 
1 NTIYPER= TYPE (0 OR LOR 2 OR 4 OR 5) OF DISCRETE RING NO.( 2 
1 NTYPER= TYPE 1 OR R OF DISCRETE RING NO.¢ 3 
0. 3000000E+08 S$ E = YOUNG’S MODULUS OF RING( 1 
7. 125000 A = CROSS SECTION AREA OF RING 12 
7.535500 Tey = MOMENT OF INERTIA ABOUT Y-AXI§ (SEE FIG. ON P.70)( 1) 
2. 375000 Ix = MOMENT OF INERTIA ABOUT X-AXIS( 1) 
0 IXY = PROBUCT OF INERTIA( l 
-1. 781250 El = RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)( 1) 
0 E2 = AXIAL COMPONENT OF EN eta en. i 
.7501E+08 GJ = TORSIONAL RIGIDITY 
0.7330000E-03 § RM = RING MATERIAL DENS TY G. ALUMINUM=. 0002535)( 1) 
0. 3000000E+08 $§ E = YOUNG'S MODULUS OF RIN 
8.250000 A = CROSS SECTION AREA RING G( 2) 
11. 69800 Ty = MOMENT OF INERTIA ABOUR IS (SEE FIG. ON P.70)( 2) 
2. 750000 Ix = MOMENT OF INERTIA ABQUT Ko ANT S( De) 
0 IXy = PRODUCT OF INERTIA 
-2. 062500 El = RADIAL COMPONENT O £4no ECCENTRICITY (SEE P. 7@ geen 
0 E2 = AXTAL COMPONENT OF RING ECCENTRICITY ( ys 
0.9249000E+08 § GJ = TORSIONAL RIGIDI Tx¢ 2) 
0. 7330000E-03 §$ RM = RING MATERIAL DENSITY CE. G. ALUMINUM=.0002535)({ 2) 
0. 3000000E+08 S$ E = YOUNG'S MODULUS OF RIN GC 3 
8.250000 A = CROSS SECTION AREA OF RING( 3) 
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wy = MOMENT OF INERTIA ABOUT Y-AXIS ( SEE PiGaON P.70)( 3) 
IX = MOMENT OF INERTIA ABOUT X-AXIS( 3) 
IXY = PRODUCT OF INE (3 
El = RADIAL COMPONENT OF RING ECCENTRICITY ( SEE P. 70)( 3) 
E2 = AXIAL COMPONENT O R NG ECCENTRICITY( 
eg = TORSIONAL RIGIDITY 
RM = RING MATERIAL DENS 3° bree ALUMINUMS 0902535) 3) 
K=ELASTIC FOUNDATION MODULUS al G. EB B/IN* 3 YIN TH 
LINTYP= INDICATOR (0, 1, 2 FOR TYPE OF LINE LOADS 
NLTYPE=CONTROL (0,1,2,3) F FOR" ry E OF SURFACE LOADING 
NWALL= INDEX 5 gmen 5, 6 al R WALL CONSTRUCTION 
E = YOUN S MODULUS FOR SKIN 
= POISSON'S RATIO FOR SKIN 

SM =MASS DENS LLY OF SKIN CE. ALUM. 90025 LB-SEC**2 / IN**4 ) 

PHA = COEFFICIENT OF THERMAL EXPANS 
ANRS = CONTROL (0 OR 1) FOR ADDITION OF SMEARED STIFFENERS 
SUR = CONTROL FOR THICKNESS INPUT (0 OR 1 OR -1L) 
ARE THERE STRINGERS (PLEASE ANSWER YOR N)? 
ARE THERE RINGS Cel SE ANSWER Y OR N)? 
K2 =CONTROL (0 OR 1) FOR INTERNAL OR EXTERNAL RINGS 
E2 = RING MODULUS 
U2 = RING POISSON RATIO 
RGMD= RING MASS DENSITY 
IS THE RING CROSS SECTION CONSTANT IN THIS SEGMENT? 
IS THE RING CROSS SECTION RECTANGULAR N)2 
DZ = ARC LENGTH BETWEEN ADJACENT RINGS ( CONSTANT ) 
T2 = THICKNESS OF RING (CONSTANT) 
H2 = HEIGHT OF RING (CO NSTANT) 
DO YOU WANT TO PRINT OUT TH CC J) AT MERIDIONAL STATIONS? 
DO YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 
GLOBAL DATA. 

ST = PLOT OPTIONS ufos -1=NONE, O=GEQMETRY Le vaw 
OB = STARTING NUMBE CIRG. WAVES (BUCKLING ANALYSIS 
NMINB = MINIMUM BROF CIRC) WAVES” (BUCKLING ANALYSTS 

B = MAXIMUM NUMBER OF CIRC. WAVES (BUCKLING ANALYSIS 
INCRB = INCREMENT IN NUMBER OF CIRC. WAVES (BUCKLING) 

EC = NUMBER OF EIGENVALUES FOR EACH WAVE NUMBER 


E = PRESSURE OR SURFACE TRACTION MULTIPLIER 

TENE = TEMPERATURE RISE MULTIPLIEK 

OMEGA = ANGULAR VEL. ABOUT eae OF REVOLUTION a 
NUMBER OF POLES (PLACES WHERE R=0) IN — 


CONSTRAINT CONDITIONS FOR SEGMENT NO. 1 1 
IPOLE NODAL POINT NUMBER OF POLE, IPO LEC 1) 

AT HOW MANY STATIONS IS THIS SEGMENT CONSTRAINED TO GROUND? 
IS THIS SEGMENT JOINED TO ANY LOWER-NUMBERED SEGMENTS? 
NUMBER OF POLES (PLACES WHERE R=0) IN SEGMENT( 2) 


CONSTRAINT CONDITIONS FOR SEGMENT _NO. Z 2 Z 2 

AT HOW MANY STATIONS IS io SEGMENT SO Ree O nD TO GROUN 
INODE_ = NODAL POINT NUMBER CONSTRAINT TO "OR 2)" INODE 1) 
IUSTAR=AXIAL DISPLACEMENT CONSTRAINT (0 OR 


D i 
ICHI= MERTDIONAL ROTATION (0= =C IN ED 2 iy OSL 
D1 = RADIAL COMPONENT OF OFFSET OF GROUND SUPPORT 
D2 = AXIAL COMPONENT | OF OFFSET OF GROUND SUPPORT 
IS THIS CONSTRAINT THE SAME FOR BOTH PREBUCKLING AND BUCKLING? 
IS THIS SEGMENT JOINED 26 ANY LOWER-NUMBERED SEGMENTS? 
AT HOW MAY STATIONS IS THIS SEGMENT JOINED TO PREVIOUS SEGS. 
INODE = NODE IN CURRENT SEGMENT (ISEG) OF JUNCTION, INODE( 1) 
JSEG = SEGMENT NO. OF PREVIOUS SEGMENT INVOLVED IN JUNCTION 
INODE = NODE IN PREVIOUS SEGHNT | JSEG OF JUNCTION 
IUSTAR= AXIAL DISPLACEMENT LAVED. L=SLAVED) 
IVSTAR= CIRCUMFERENTIAL DIS AN EMERY O=NOT SLAVED, 1=SLAVED) 
IWSTAR= RADIAL DISPLACEMENT O=NOT SLAVED, 1l=SLAVED 
TCHI = MERIDIONAL ROTATION (O=NOT SLAVED, l=SLAVED 
Di = RADIAL COMPONENT OF UNCTURE GAP 
D = AXIAL COMPONENT OF JUNCTURE 


GAP 
IS THIS CONSTRAINT THE SAME FOR BOTH PREBUCKLING AND BUCKLING? 
DO YOU WANT TO LIST OUTPUT FOR SEGMENT( 1) 


Aas 


DO YOU WANT TO LIST OUTPUT FOR SEGMENT aa 
DO YOU WANT TO LIST PREBUCKLING RESULT RING FORCES? 
DO YOU WANT TO LIST FORCES IN THE DISCRETE RINGS. TE @ANY ¢ 


oo 


2. ACESNID INPUT DATA FOR ALL CASES 
The input data for ACESNID must encompass all waves N 
used for the BOSOR4 model. The waves are specified by 


NSTART-NFINISH. 





VIRTUAL MASS FOR RING STIFFENED CYLINDER WITH FLAT ENDS, N=0,1,2,3 
0 3 1 1 0 _/NSTART,NFINISH, NFREQ,NVMASS .NCHECK 
12 6] lee Leu NCYL NRITE NWBOSG ,.NSYMF 
2 3-0 {N RDER NFEN DS yN CENT .N FOMET. NOMIT 
9.59684E-5 5.83354 0 5”) RH SOUND. ERR 
: i Like 1} Oil 1° 1 OUTP i OFL.9 
de: ‘ 9 )SETS JSGBEG,JPTBEG,JSGEND,JPTEND LEFT 
eer [SEES JSGBEG , JPTBEG ,J SGEND JPTEND CENTRAL 
1 1 1 9 /SOURCES JSGBEG TP EBEG SGEND ND LEFT 
2 1 2 27° /SOURCES JSGBEG, SPTBEC, SSGEND> SP TEND CYLINDER 
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oe PIFLASH INPUT DATA FOR ALL CASES 
BOSSY100-103 are the shell mode files produced from the 


BOSOR4 runs. Note the deleted shell modes that were purely 


meoemsional for N = 0. 
4 00 11 /NUMBER NTORSN NPTM NSYM NSYMP 
m2 2 | (NWETSEG K),NUSESEG(K) ,K=1,2 
386. /GRAVITY 
BOSSYLOO /SMF,N=0 
BOSSYLO1L /SMF,N=1 
BOSSY102 /SMF;N=2 
BOSSY103 /SMF:>N=3 
22 30 30 3 f (N USE(J) Jal NITEMS ) 
e233 { KORGS K) = NKORGS | 
m46 7 13:16 15 16 17 18 19 20 21.22 23 24 25 26 27 28 29 30 
meg) 4 5.67.8 9 10 11 12 13 14 US 16 17°18 19820921 22 23 24 25 
DG 27 28 29 30 
12345 6789101112 13 14 15 16 17 18 19 20:21 22 23 24 25 
26 27 28 29 30 
12345678 910111213 14 15 16 17 18 19 20 21 22 23 24 25 
26 27 28 29 30 / (JUSE(J,K), J=NJUSE(XK) 
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CASE 


Sa Daa 


FOR D 


SAP IV INPUT DATA 


4. 


Vv 
4 
Sy 
rs) 
O 
= 
M4 
aa 
Y) 
oH 
GU 
¢ 
hy 
Vv 
yp 
cS 
oc 
Vv 
<a 
ou) 
kw 
Vv 
Ew 
© 
4 
Vv 
Cc 
Vv 
mM 
Vv 
re 
om 
U4 
G 
pV 
© 
wT 
m 
Cc 
A 
> 
O 
om 
ro 
O 
Wy 
Vv 
Ae 
e+ 


Note that only half the doubly overhung beam is 


model. 


modeled as indicated by the beam lengths at the node, i.e 


- 44.81250. 


node 21 
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Dis PICRUST INPUT DATA FOR D = 4.2 IN. CASE 

The circumferential location of the substructure 
attachment point is specified by ANDEG which is the angle 
ccw from the shell global Z-axis. The meridional location 
is specified by the node 21 at frame 18. The angle DEGROT 
specifies the orientation of the substructure local z-axis 


to the shell Z-axis. This 1s indicated by a ime PIgure 5. 


mets Ort 1 oro oO 

pete 1! 0.10.) (QUIPUT FLAGS 1-20 

5 0 0 /NJUSE,NHWSOB,NHWBAR 

i271 tO. O [NIPSUB , LBOSEG , LBOSPT , ANDEG 
180.0 / DEGRO 

i? 30405 (Cs 


Wsh( J )ji,J=1,NIUSE) 


8 9 


6. USLOB INPUT DATA 

In Phase III of ELSHOK the USLOB code performs the 
integration of the incident pressures using a Runge-Kutta 
scheme. NCHARGE specifies which type of UNDEX is to be 
represented. For the tapered charge nine discrete pressure- 
time points are used as input. The numbers following the 
charge weight of 352 lb. are the constants used in equations 
(1) and (2) that empirically describe the explosive HBX-1l 
(order-Kl, K2, K3, K4). Note the 0.0 for surface cutoff 
effects. For the empty shells NSUBS must be set to 0. The 
integration is performed to compute velocities at all the 
nodes that make'up the shell and substructure but only.a 
couple are requested as printed output in NPTSH and NPTSUB 
as indicators that the program ran successfully. Nodes 1/01 
and 2/21 are checked on the shell and node 1 and 21 are 
checked on the substructure. The value NQUAN = 3 requests 


the z direction velocities as Gupene. 


TAPERED CHARGE 
1 /NTIME .NSKIP ,NCHARGE, NQUAD, NFINE,KOUPLE, NSUBS 
3° dp-05 123. 0625/856. 8123 0.0 /DELT, xHGAD- wate SURCUT 
NS 
0.2 0 .00053 350.0 90243 305. .00294 160.0 
Bo 0as O18, 0298822960" 988820728: 9 88242 895 Bx80RBHE BELA? xveurs 
Hee io. | O1 9 1 
101031061010 /OUTPUT FLAGS 1-20 
p (NETSHL 
1 1 2 21 /LBOSEG,LBOSPT 
a (NETSUB 
1 3 21 3 /NSTAT,NQUAN 
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7.  PUSLOB INPUT DATA 

This data file generates the velocity profiles that are 
plotted on Tektronix devices. The time step DELT and the 
interval of integration specified by NTIME and NSKIP must be 
the same as for USLOB. This input requests velocity 
profiles be expressed in terms of the shell global 
coordinate system. Shell output is requested by segment 
node and angle. Each substructure point is specified by 
node, number of degrees of freedom at the node for which 
velocity profiles are to be produced and which degree of 
freedom it is produced for. The tip is node 4, one degree 


of freedom in the local z direction. 


1/NTIM.NSKIP ,NSUBS , NTEK, NCARD 

a eaecoe +665 120/DELT , XMULT , YMULT 
VELOCTTY PROFILE D=4.2/TITLE OF PLOT 
TIME tis “InN /SEC* 
16/NPTSHL-W 
Venu 1/04 
Taos 
t aed: 9 180.0 
} Ard 3/939 280-¢ 
2 3 1 2'0.0 180.0 
VeaTH 2/07 
Ghat 3,999 180-° 
2 12°12, 9.0 180.0 
2 tg 129.0 180.0 
2 21 1,5.9.0 180.0 
327 1,2,9.0 180. 0 
y- ATH +1 80.0 270.0 
V_Annae? (Od 
9230) 1 90.0 
Vontne? 1a 
27 11 90.0 
V-ATH 2/12 
20 ee 9050 
7 Aa 2/16 
Jia el 690.0 
yap eies0 0 

= 3 /SUB_TIILE OF (PLOT 
5 aTH eels /LBOSEG , LBOSPT , NSPHC, NANG , ANGDEG 
16 (NPTSHL-U 
V-FWD 1/01 
t twh t/Q4e 
t-Bup 4/830 180-2 

1 9 1 20000) 18000 
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APPENDIX B 


PROGRAMS USED TO PREPARE AND PLOT VELOCITY PROFILES 


PROGRAM PLOTCONV 


THE PURPOSE OF THIS PROGRAM IS TO READ THE PUNCH-CARD 
FILE FROM PUSLOB AND CONVERT_IT TO A FORMAT WHICH CAN 
BE ACCEPTED BY EASYPLOT OR DISSPLA ON IBM 3033. 


QAAAAY 


cane ay REt eS 32), VELO(166,2), LABEL( 32), TAG(32), LIITEE@aze 
RACTER*40_ LTITLE 


C 
sci 

NCURV=32 

DO 10 K=1,N 

READ iy -9995- ae »NRE S ,DTRECS 
READ( 11,998 G(k} 

READ( 11.997 LABEL yt RECS) 
CONTINUE 

FORMAT( A40,15,1P ee) 

FORMAT( A10.F10. 

FORMAT( 1P7E11. 4 

DO’. 20 L=1,NCURV 

T=0.0 

LL=L+20 

DO 30 M=1,NRECS 

VELO(M,1)=T. _. 

DTRECS1=1000. 0*DTRECS 

T=T+DTRECS1 

VELO(M.2)=VEL(M,L)/12.0 

CONTINUE 

DO 40 KK=1,NRECS 

WRITE( LL, 996 ) VELO( KK,JJ),JJ=1,2) 
WRITE LL, 995 LTITLE( }) NREGS, DTRECS 
WRITE(LL,994) LABEL(L) ,TAG 

CONTINUE 

FORMAT( 2F15.5) 

FORMAT( A40,15,1PE11.4) 
FORMAT(A10,F10. 5) 

STOP 

END 


94 





R 
10 


20 


oe” REQUIRED TO PLOT 
2( 166) ,VEL2( 166) 


0 3 8. Ne 
6. 9 5. ) 
TIME (MSEC)S$',100), 
ViuOStEte FAGBHE 6°80, 0. 2 


Y 
2731 ORT. ve Ahi 


10 
: Sueeie! 
ED SHELLS 


te 
Se 
ester Sel lele]etule la 218 Late oO 


ee, CHC Whe! 


S 


CS Fico eles 7) 


BBY 


9", 100) 


rcs 
Scaes > OOO ORMMHH 
HARM WMnCaACrsSsaasaswM 


aaa aVEL, 166,0) 
TIM2, VEL2,166,0) 


ce 246.0445) 


GR 

RES 
LE 
EN 
DO 


Soaun 
CIFIC 


ONS 


APPENDIX C 
SUPPORT BRACKETS 


PROGRAM USED TO ESTIMATE AXIAL FORCES IN THE 
PROGRAM FORCES 


-_~ 
Oo 
oO 
=~ 
— 
oa) 
fx} fa 
At 
A 


BB 


9 
9 


on) 


) 


ORT BRACKETS 
166 
Jy 


TIONS 
2 
1 


PP 
E EC 
5 VB 


SU 
FL 
16 
0, 


WAM 


hiv 
HA 


Ti 

Tat 
“4 XL 
*4 XFO 


OOO 


CONVENTIONAL CHARGE 


LBF 


1,2) 


MSEC 
(2X0 LN) NN 


) *EMOD*AREA 


, 


XL 


ON 
ONIAIQIQNQUQHE & -o~ 4 


HOMA Mt St AW HA 
HOMAIM + + FARE 


P 
166 
EMO 
RCT 
ai 
iis 
TAPERED CHARGE 


i 


ESTIMATES FORCES DEVELOPE 
9 
t 
4 


OND HHH HZ0OR 
fret IO O MMOMOMOHZ 
ASSOSAzOUMN 





} 
UETVNO TFT NO DWOMOWMAUNE.cOr.wWr10 
PN SAS NNMNOM AAR Et ART InNw! 
wad eC ce! wu. cu ge eee meee once OPCs ees 
PANS OMNDDWAADAID CM oun a 


Ber ATOM MOMANMMAKT I Nowe to: 
MOAN TMMAKHMFHRMRSOWMMINANC: 
AtOONAWNTOTNIM AM NIN AT: 


PMN AOM OM ~SF UO NOT OFAN sei 
teen an te 


{x) 
HOOOCOOCCOCOCCO0C0C000008 


ONDONTODOORANMMNOMDOHAGM 
2IDNIENFUAOPS ON wk et klk ltl tlt tll 
Wien © 6) 6 eens Sse stag 
TSNoOooco0co0o 


NNO DAA OUMNAROtNAORo- 
FANNOM SFA NAAT OS aWOTasNS 
CIOS 6 ew MONA 6 1 ON fe ete 
© COrmd . « « OCOOMD wOO . 


fu -rANnTFOONDOCAOWRTOROWNG 
PONS ASTIN AAR A Anmnomw Omanwn 
ATONKDNOMMANMALSTOMawMsteuneae|g) 


ORT AN SIT BORA MIM NNCOUNa a 
teeee ve 


N 
>] 

¢ 
{x} 
HOOCOCOCOOOCOOOOCOC0000oO 
HOOCOOOOOCOCOC0O0O00oo000 
HQOOCSCSCOOOOOOOo oe OOeooo0 


ONDONTODOONMNMMNOMODANM 
IMACS UOMO «ww ew ew we th tht tw 


2G 


IANA CRT AWMMNM DDO MMO ANW OT OUI OW MT AMNMMNDNIDAANODWONAIMo~TOCFOO 
POD « COAQRTHOIMNIMUONANMOMOOGHT MTom ONDNO-T OTA Gsanmno~Thcoomomrm-~T 
ee owe oaqN ONO ON} ew ew tl lel lw lw oN Uw ct Cw wl ltl Ct ON On tw tl et eNO 6 we 0 le oe etme ee s 6 e 
LAAN UI) AON OOM ODOMO MOMDOoT MmNOUMO . TONROWM MRANTNROO 
MODOOMOOMELTNNAOMADOMMNAYYE® ANWDWANMMOAMWNDANROAOOMAAONNOAM 
ISPS WICIMOWAIURT WOON 0010-F OM ND 00-F OD NNMONANT MR MWR TANANNOT OOM Nou 
MANOFAOTIN TOWN WOAMANTFIOANMMHODMUODANONMHTINOMMATAOMONON 
PAN WIN ARO PR RIRIM DOVER QIMET ANGI AIIM AN MM ADNANGBT MPR WIR TF AM WNIOIM MAN AST NAN 
'eoe 8 Beeeet + 8 ' ee ’ 


DDD D COO DCOO COO DVDOCOO DOO ODO COOOOOSOSOSOOOOGCOTOOO OOOO OO0OO0O00 
DDD DO CODCOD OC OO GOOD OQOOQ QO QO QQ OOOO SOOO OSOAADAAARAAARAAHAAAAAAA 
OOOO DCD OOO ODO OO OGD O OOOO OOO SCO OOO SOQOHAADDTDRAAAHAARAAAAAAAHAC 
DOD DCW COO O DO OO CO DO OGVOOSOQODOO SOOSOSOQOOOAADNAARAARANAAAAAAAAG 
WWDONTWDONATODO GET ODON~TWODOOKT ODONTODDHAMNEHAAIMMNEHAAMUIES art 
PUES OMO ASR TOM NO ANMNORONDANEGM-TUIP WHO AGT NOM DOANM-~TOMANON 


QIQIQAI(IQNQUmmmn MONON MNO FTP Ft TT TT I 9 U9 90910. 00 00 00 0011000 COPA PRBS PRPS PS PRP 00 


NDANOTODOWNIMANNADAAMMONAMNMHMNOMDOMAADNNAMmeTODUWM NODAMMANUNNO 
WOMDWO MRR NANT ANDI AMO AR ARTF DM Or AT TRHMUOMEMHOMDO ONAMAAT MH 
© CO COCDCIMDMAWMd 2. AMN WKN ogurtRewoomnd . TITNOMAOAN .. Ono won nw 
MwOdcomMnN . oN « «CNIS 6 5 oO .0O we, SCAU) oe oc es ee GN .coOUnNUD . UD LON ew 
wWOor-wouw oor Slips De ead _LODUIR MAROON OnsS ODOUR RWOM XO 
ADDM APS ODDO MNAN AUR UMD OWUNOMOOMNONODOTOMAMNNOAAAMMMOA 410 
DADNDO MOT FTHOMM HADOMANWDOTFODHAANAT OUINOO-FTFDAAKHAMTFOOWMO OTH 
et Ses AD ART MOON Se ATO AON TNS HOO OMI NM At eet Ot SAN NON oe 
166 886 86 6 € 8 8 8 & & 86 ’ 8 ’ 8 to 6oe6 8b 6 € 8 8 86 


ODDO D DOOD OOO DODO OOO GTO OOO COS SOOOOSOSOOSOSOOOOOOOOOOO0O0900 
DOOD DC OC OCQOQOO0OO SO COQ QOOOOOO OOO OOOO OQVOWDAADAHAAAAARAARAAAAAAADS 
DODODDOD DODO OO GD OOO OO QO SCO OOOO O SOOO GDOADAATDARAAAAAARAHAAHHAC 
DOC O DCO ODD OOO OO QO QQOO QOOQOO OS OO OQQOAAAAARARAHAAAARAAAAAAHACO 
WDMOaTWWOAATODON~T ODO ORTODONTWODONTWODNAMYVIEHAAMMNMSAAMUIS oar 
TUES ODO AM TMDOMHNHOANMDNOMANHANUMTYVDHAOMAALTMNIWOM DOMAICNTOP.DNHON 


RIQIQAAAQ MMMM MONON MNO STF FP TTP MUU U9 0109 0 00 00 0.0 00 1.0100 LOPS PS PRBS PS PRPS PSP 0000 


MAMNOMMO WOMAN 
an OL SOS STO) 


NoMOFOSORRONON 
IMMINMNNON SINGS 
6 0 $6 


ejalealelelalelelelelalale) 


ODD DOO DONHAAAARAG 


CHINHNONTOORMNOSN 
DMWOMODOMATANMMFIN 
ie tel a 
Nope ae oemos 
ArT WOFIMMNMAMOOnm+t . 
WANOAMNMNARHM~t 
St ANT COD COP.~TFONACO 
’ ’ 


OOOOOCOOOO0000 


OOOOOOOOO0O000 
FTowodesTwowodn-Twoo 
MTS OHO dens w0or 


COCO DO CODD DO HAARAAAG 


oF 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22304-6145 


Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93943-5002 


Profeseor Y. S. Shin, Code Gesa 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5000 


Professor R. E. Newton, Code 69Ne 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5000 


Department Chairman, Code 69 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5000 


Dr... Neen Cs al 

Defense Nuclear Agency 

SPSS 

Washington, D.C. 20305-1000 


Dr. Hewesevin 
Defense Nuclear Agency 
Washington, D.C. 20305-1000 


Dr. Meee saron 
Weidlinger Associates 
333 Seventh Avenue 
New York, NY 10001 


Dr. Ranlet 

Weidlinger Associates 
333 Seventh Avenue 
New York, NY l1OOOL 


98 


No. Copies 





EO 


ea. 


ivan 


3 . 


Dr. Andrew P. Misovec 
5202 W. Military Hwy 
Chesapeake, Virginia 23321] 


Pans Daeumang, Code 1750.2 

Hull Group Head, Submarine Protection Div. 
David Taylor Naval Ship Research 

and Development Center 

Bethesda, Maryland 20084 


Dr. Huang, Code Rl]4 

Naval Surface Weapon Center 
White Oaks 

Silver Spring, Maryland 20910 


LT G. F. DeConto 


47 Water Street 
Sandwich, Massachusetts 02563 


goo 


[09 \A) 





















































and 2 t . 
¥ | a So 8 Fi ft a er er er ee i . r ‘ ' 
al i ad a | y a edory | ik URL ae ba oe a ae "cee er ; 
SF oibivaplecé ian ee one ‘ Ce el eT 2 r a” « K ' n r p 
MEiit ot ot ee kt 8 A Pte Ee CNS "fe ayy ATS | ter eo Fe a I ts y 
Le eek oe ee ee Le ee ae hatte wig, aay er ne ee Oo Cn ¢ a L - ] Y h D234435 
Os chi teble ty +7 aL 448.1 " r ia: é eC ae 2 t es 
eet Ra el Or ie eat oY rin ye Se a ee P 48 ah - 
oF hn re | oud x iP cer 45 ' co i r ei OO uidy a ' ‘ 
tar LTA LT aE Me Ha Ee A parametric study of elas S 
al ek i i ert sa! Cho ee ee ae a ee Ce 28 ee he t ‘ r v ponse oO 
re Ca ar ee | Or ae | a) re | [an | a 5 F P 
2 tad {fa ree Sree? ee ae ae. 1; Fi ee ee | ee oe tf 
t 1 Te ae et? "1 a Ay ij ' | , F 4 ¢ S f 4 3 
iy Rie alr 7 ‘ Fry rad i ae me he la? ¢au ar i ‘ 
DATO ee Te Set ot ee Me ay LT ae Me Le Ly Be A i | ai Y | Pi 
eT Lae a n Cy at oe. Te oe ee Ed rue i ee f . 4 
fe ee £ Oe D ‘ 4 ee Oe ee eT | aay) re ‘ 0 ti A ‘ ss ' i 
Pie Tid Pee ’ oO RT 7 ‘ ar ay ptt L Fi o ere, ' ts Lae a 
7 ore ar hed cy ct an ee ee he | rl re a H 
n erode ard oe ee del ne ‘ ae; 5 
f Oe Ue ee ae a an) 7 4 ‘ 
6.6 4 Orff Cee ar Pref AH 8 ad 
ae fn a ee | as } F i ] ‘ /68 O 
ee ee oe ep . F it 
t eT es i et eee en oe on | 5 i are | iil 
Ole ee FP ote 8 aT s ¢ " 
. ih Pe ie 3 ry é {~ Sn | My a) id a DLEY KNOX LIBRARY 
£ a ee bn ae ey fe | fy c ra a t ' ‘ f 
Bi 2 aha ! j . ia f ‘ r , 
FoF Fi Ca) re) 4 r ‘ ‘ ‘ iar) " 
, (an é | r ‘ ne 
$ of r ‘ S08 ' i r t t L L 
Tee ins" 6 i: A i i 2m . i o , A 
¢ s 4% 1 Fi a an ‘ ‘ ‘ fy 
- P ae an a a Fi o 5 D : Fl A 
P) # ae é e ed uy ‘ A 
CY (a te ed ‘ fT r f ] LI ' + 1 
fay) PF . *. ro il ' eT bd Uy ‘ a) 
Py at ‘ eS f a ¥ aay ‘ 4 
U im i L ni Pa Ce LJ s ' : 
fi an ft of rd ' r | ie | e 
fF toy Ce ‘ ‘ a s ‘ a. 
r F Pee ee | D ‘ a r D ¢ 
ae i ga ‘oc uc f - r ery » { " ' 
r r ‘ “ ‘ ri Py ra t t ae 1 
4 rae D <t rd Cars e o n 
A ¢ Fey Feed ima © 5 r f 
au ‘ a ‘ P Fa at 
* Pi te fi * Cd oe @ wo. Ei r oi 1 
p Fi a i , r] r ra ‘ r 
P ¥ a oy 3 ae) oes o fF i f 
v4 y r “d iifae a? ers. ' a Tr ‘ r ) 
‘ P PY 2 r r £ ‘ 
Fe a a 7 Py » Ce. | Sf Yan | D ‘¢ 4 ‘ 
f # yd , ore J - a ¢ i 
* Lr) a 4 * A 1a 2 ‘ a P 
Q 48 4 a] v. ® 
‘weet, tl Ui ry o Ce ‘ i . Z ‘ 
' dl ad a < ‘ ee rT J] i a ’ . A 
r Pr Cer * £ Pe 7 
aie ’ r f ‘ f F 5 
ro * é & ra ’ ba U . i] rs 
a Py te . - : 
] Cs a. 4 
e . J ’ L 
. Pa * Py Py fy 
Py | 
‘ é 7 4 ae) oa 
ry i an 4 ‘ é 
4 a. f ag ‘ 
iT ' * 
4 ¢ ti ‘ 
7 d ’ d 
fs a LJ 
5 n A f r 
f , x 
ro Lf ‘ 4 
‘ 
n Ps 7 * I et 
os * r 
Ps 
n A 
r LJ ~ 
. r 
‘ 5 
rd "4 Y 
a * 
‘ 
P ‘ 
, 
: P 
a sa 
® 
f 
: . f 
F - 
* a 
. 
J ’ 
LU Li 
a Ss 
Fd 
f 
5 
n 
A 
J 
- 
A 
r 
, 
2 5 
‘ 
. 
5 
5 
y F 
rs 1 A 
eo. 
oT P , 
Ch ‘ fe A : 
5 
. my or F 
o 5 
A * 
7 ‘ C i 
. e 
7 
“ a A ‘ 7 
. 7 
% - F 
7 s 
/ ‘ 
T i 
slid . * as " 
7 bets Ms 
p . ‘ 5 F = 
b 
2 . 
c A , 
F e -« , : 
n , , 
a . ae rs 
P F . * o 
“is : 7 
y F ak 
Po ry Cra ° P e 4 : 
* ale D i 
. ~ na ‘ ® a 
. Ps ® n 
sd L . 
ee * a . ba L ba be s cf ta bd 
x rT : P 
. iT e iy + ; 
. Py e ie e 1 s e P P - - Ps 7 
sn ° ad Li Sue® a el % 4 
re4 e Ci % a as 4 + = Ps . 
Py C oe Ls | itm Py : i * 
a L hd ary & * ce s oy r e i LW 
b i? M, e ‘ wrath ir = if t a’ wh Tel ‘ ' A . 
‘ L Cf . A B bs ? Lm “a : . 
Fe Vy ht ek , a i o rj 4 A 
A bed he kel ™ “<i rN 4 F , e f 
al" 2 ot L a 4 a She] Ce ee cr Py ry 
eh tl Ba | Be ee td tit hed al py r ‘ 7 . 
Od hal Le a oat Meek ol LY b i =% ' 5 
¥ 4 Tals ie id i ae eo fi a a S r 5 " 
a bit il ak % r La ea eS ‘i . r 1 . 
& ¥ ab . La | oe Sh CLE e = an] at ' 1 
bhadals toltien beat ak eh 4 t yu: n “8 aa r D 5 
ed 5 da Oy * Ki oo. Tat | r 
a Ti ek, hag Oe hl F P ' ae . Py Ser tr o fo A . , 
; ed a "Et Mery L Pt A on 
Ct oh bait Toads Td P n - ry : n 
ty > i & be 4 * s ee ir 
| av T . te! § a | wed, , ] . 
he er & ea | ah ‘ ‘ o 8 r ‘ 
ide Le he he it el bel eh, a he A Sea Tot ie | r i 7 
a Lave a or a afete ar Py O / 
sewderh sp* er ee ry Pl p ae : ry 7 
Oe | eb | ti wey a | * ‘. o & S OI Y : 7 ‘ 
oe a Mal | “. dale Siete Pa Lot ie Sel” t nai ied n ‘ *ua. ed r . " , , n S e 
wegagd iar ee i ee ie rs . Le rt Yan be a i so 6 ry . 
t r Ps ' . te A f 4 ry ri \ te | or] ae eyes ar | Pa a ‘ : 
+3 U My ' . ih ee ee | KMhageg ee ee er cd “lie H iia ‘ Ps 
ta Des he halle tel eh held eo i Vigce ng raPl Portion, Cae oa ae n Px O ' . : 
' + a“ ha ee tel ae oe | eh htedk i ee ey ry a. PT P P F n - t 1 . ‘ . 
Wiel | tg ry, Why, ’ ht ob Le ee | rm i i i rat Par) A F to. = P 1 . 
UE Ot ek der | uy he eed Po a PL ee LA] « Fr 4 r eta 8 ‘ r Py n 
tol Tedidied ol UM iol i aa b Le rant r) P : , ‘ A 
Li Ln ae: ¥ en Fi ioe ‘ tear te rt Py 7 
use Hy iY \ pees uJ bt dey kat al i : : . 
ih Tedd thd bd re ‘ bus ek a ad 3 , 
“yn at * a oes a r rl : ' aa 
Ce ae roe c Li an e i a . ry ' e fo 
rs ea te Te Lr tel ge at an ' Py D f : : . 
Weary 4 Ww oo a eee Lr a a | F n 
ee ee oe % a hi LU "vu 
lh a eed ted To ek s P we ebit x "t % t. r » ¢ ui rT ® p 
FIP oe peeve tN 4 Lf Y Te a OTT feos Fo | gbe A ’ a 5 , ‘ ' ee : 
EW Ayebres oy yobs ure ye Uoraderyseclsyet Oh eT ehigs 5 | A ' A A 
ee} a eh Oa oh a UE Th ey an arte er i fT ro ¥ " " 5 A F ‘ ‘ 
ae il | eal ih) ile far ra : a) ¥ es ee a | a ea oe ® 
Pert hee ed On bt | fet Uw g 4 eo eS dt % 1 PT ir P . 1 
Vege de ee ee it a Oe oe er ee Yi 4a re 3h) Riding 
th aga Weetate bs wel. A'@ Viued Pa oun ‘wg LG " P ' 
‘ i <4 ee ee Od) ee ee oe ee a Ce hl ee ee eu va tr 8 + f 
i a tie ee tee a, eee kT TE it es t tpt eh et i n 
arafury ey Be -qmenta gt bob ee 2 er i has | a tet a | ait ar toe a 
OTe Sey a PP beviven iy irs EAT) rele "i a%’¢ sat 1 p r “Pa ; , 
ee i ee te aT fe eae 2 oe f ' ae Pe ae ek an] 4 : p . 
di Del tad ge % rt & @D4aPget ig i Pl ar ‘a.m ¥ neh } D 
CU ae er Vibrks © OF aHTY Bees 8a at Pr a ee rae ear rl y ; 
Oia utere ry a Loe on ee Che a 1 te oe i i he ee t‘y 6 U uJ 2 
SMe ed a C et. ee ee ra a te ar a ee Per F ' : 
Oe eae wea yay Uryrtih ay Ca an er Pe i ee ’ f 
ee ae ge Pa tt €wlethun.g g wie % we bs hay an! A n F 5 ml ’ ‘ 
a a) a ee 2) ae © ee | vy OW ba gg r lwe } : a) ' H , 7 
FA | nr 5 \ Ut Hy Usten es a) a Oe 1 r ay hey rae \ ar | o 
ae 4 oo on Tr | i fie a4 by es i "Ci iar tar i 1 
- a es eeu ay Ci i Pan ot ee P ' . - 
1 is U4 L 5 » i i r Genwy,t re p ’ A 
Sy et BN aCRsE 4 L B.rURt on a: — UTS ‘vy a 7 ‘ TT t ‘ ‘ . 
wi Lite « & yl L ~ Oar a a ro n 
Ue acd | eo VP OP sy A i ¥ rt (Sn ree .. rl A 4 5 Fi : 
t auogre’s %¢ “td P a 9 Mi ry a ) i. L : t 
if Pa Wor or rl 5 e'hoe 3 oa Tee ] Py 
L ft ion ie ay s 





